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Abstract 

I provide a straightforward proof that a simple harmonic oscillator perturbed by an (almost) 
arbitrary positive interaction has a perturbative expansion for any finite-time Euclidian transition 
amplitude which obeys the following result: the difference of the sum of the first N terms of the 
series and the exact result is bounded in absolute value by the next term in the series. The 
proof makes no assumptions about either the strength of the interactions or the convergence of the 
perturbation series. I then argue that the result generalizes immediately to a much broader class of 
quantum mechanical systems, including bare perturbation theory in quantum field theory. The case 
of renormalized perturbation theory is more subtle and remains open, as does the generalization 
to energy levels and connected S-matrix elements. 
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I. INTRODUCTION 



Perturbation theory has historically been the most useful tool we possess in analyzing 
quantum mechanical systems for which we have no exact solution. In particular in elemen- 
tary particle physics its predictions have been confirmed to astonishing accuracy, for example 
in the four-loop calculations of the anomolous magnetic moment of the electron and muon. 
However from a mathematical point of view it is clear that for many interesting systems 
the perturbation series is a tool of fundamentally limited value. For models as simple as 
the quartic anharmonic oscillator, the series is known to diverge for all non-zero values of 
the coupling, and moreover well-known "non-perturbative" effects such as instantons and 
solitons depend on the coupling constant in ways that are invisible to perturbation theory. 
These issues were not understood until the 1970's, after which consensus developed around 
a standard folklore which I will review briefly in the following section. I expect that this 
folklore is basically correct, but there is a key logical step which I was unable to find in 
the literature: to really evaluate the usefulness of a given truncated perturbative expansion, 
we need to know how far away we are from the exact result. 1 Indeed without a bound 
of this sort, perturbation theory is an uncontrolled approximation and its use is logically 
questionable, especially in light of the mathematical difficulties I mentioned above. We may 
of course provisionally use its remarkable experimental success to justify it, but it would be 
more intellectually satisfying to find a mathematical explanation. 

II. THE STANDARD LORE 

The basic assumption of what could be called "naive perturbation theory" is that the 
transition amplitudes of a quantum mechanical system with Hamiltonian H + gV can be 
expanded in a convergent power series in the coupling constant g. Moreover it is not hard 
to prove that a power series in g that converges to the amplitude for g < e, with e some 
positive real number, also defines a holomorphic function in the disk \g\ < e. Thus we are 

1 In several ancient papers [l|, 0, S] it is intricately proven using functional analysis that for some fairly 
general systems the perturbation series for the energy levels and eigenstates is asymptotic in the sense 
that the error of truncating the series at fixed N vanishes like g +1 for sufficiently small g. It would 
however be preferable to have an explicit and calculable bound at finite g, especially given that for many 
experiments the perturbation is not something we can tune at will. 
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really assuming that the transition amplitudes can be continued to holomorphic functions 
of g about the origin. It is only in such circumstances that the perturbation series can 
converge, and only given such convergence that we can make the error as small as we like 
by adding more and more terms. 

In fact for the anharmonic oscillator we can see via two separate arguments that the 
perturbation theory must diverge 0, 4, 5]. The propagator for this theory has path integral 
represent ation: 2 

D(g, Xi ,x f ,T) = (x f \e- TB \xi) = j [Vx\ x x { e~ J?"*(^+^) (1) 

The perturbation series is defined by expanding the exponential of the interaction and cav- 
alierly exchanging the order of summation and integration: 



D(g, %i, Xf, T) "=" ^ D k g k 



k=0 



D ^ K -jrJ P>*E [J o dtxA ) e-#*^J (2) 
Since this section concerns standard knowledge I will simplify the notation by sketching both 
non- convergence arguments for a single integral whose behaviour turns out to be analogous 
to that of the path integral. 3 Define: 



A(g) = / dxe~y~^ x ) (3) 

J — oo 

The first argument begins by showing that A(g) is analytic in the g plane for Re[g] > 0. 
For Re[g] < the integral doesn't exist, but we can prove that it can be analytically 
continued by smoothly rotating the contour of integration opposite to the phase of g to 
ensure that gx 4 > 0. Finally we show that this analytic continuation has a branch cut 
on the negative real axis, whose discontinuity can be estimated by the method of steepest 



2 In this paper I will always assume that the path integral is explicitly discretized and is thus just a finite 
product of ordinary integrals. However since the bound that I will derive will be true for all finite values of 
the regulator, it will be true for the limiting case as long as the limit exists. In ordinary quantum mechanics 
we know that it almost always does, but in quantum field theory I will opt to keep the regulator explicit 
(at least in mind if not in notation). 

3 The steepest descent arguments for the path integral case require a discussion of instanton solutions of 
the equations of motion. The reader who wants more details can consult @, 
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descent to be proportional to e~ J ^\ . Since the function is thus not holomorphic at the origin, 
it cannot have a convergent power series. 

The second argument is to estimate the large orders of perturbation theory and show 
explicitly that the series diverges. We thus need to define the perturbative expansion: 

oo 

A(g)"="Y, A k9 k 

k=0 

(_-\\k poo 2 

A k = L F Lj_Jxx«e-± (4) 

This case is so simple that each order can be evaluated explicitly in terms of gamma func- 
tions, but a method which can be generalized to path integrals is to use the method of 
steepest descent to approximate in the limit of large k. This reveals that ~ k\ for 
large k, with subleading factors that match the result of applying Stirling's approximation 
to the exact gamma function expression. Thus A^g k will grow arbitrarily large with k, 
regardless of the smallness of g. 

In fact we can also show that the k for which A^g k is minimized also has A^g k ~ . 
This result, along with the divergence of the series, is often used to justify claims that the 
perturbation series is asymptotic and that the size of the smallest term is the size of "non- 
perturbative" effects. However to really prove these results we need to bound the error of 
summing a finite number of terms of the perturbation series. A bound of this sort for the 
simple integral is presented but not derived in [9(. I will now provide such a bound for the 
path integral. 



III. PROOF OF THE THEOREM 

We begin with the exact path integral expression from equation (Q, but generalizing the 
interaction to an arbitrary positive function /(x), 4 defining I[x] = §Tdtf(x), and suppress- 
ing unnecessary notation we have: 

D{g) = J [Dx] e~ {So+9l) (5) 

4 ,f(x) is also restricted by the requirement that functional integrals in the perturbation expansion exist. 
This will be taken into account in the statement of the theorem below 
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The error of N orders of perturbation theory is then given by: 



D(g) -Y,D k9 k = / [Vx] e"* £ izg^_ (6) 

k=0 J k=N+l 

Here expanding the exponential inside the path integral is perfectly justified, as is moving 
the finite sum of the first N terms through the path integral. Pulling out a common factor 
and relabeling the sum, we get: 

N 



(7) 



D( 9 ) - £ DJ = (~g) N+1 / [Vx] £ ^ (-91) 

k=0 J Lfc=o v ; ' 

Now the key step: the sum in brackets has an integral representation: 

E ^^^ N'-^- (8) 



o (A + iV + l)!> ' JV! ..„ 

The reader can verify this by putting both sides into Mathematica, but I will also sketch a 
proof in the appendix. Inspecting the right hand side of this identity we see that for x > 
the integrand is positive, increases monotonically throughout the range of integration, and 
is bounded by the function t N e x . As such we can easily prove that: 



1 1 

°sE(rnf+TF<-*>**w+i)i (9) 

Now going back to our expression for the error in equation (j7j), we see that the integrand of 
the path integral is positive, and given our assumption of positive f(x) it can be bounded 
from above by replacing the sum in the brackets by , N+1 y ■ We have thus proven: 



\D(g) - g/Wl < j [Vx] I^V*° = \D N+1 g N ^\ (10) 

This ends the proof for the simple oscillator system, but the argument depended very 
little on the forms of f(x) and Sq. For example if we had considered a system of many 
oscillators coupled by a positive interaction, the steps would have been identical. In fact 
the only requirement on So is that it falls off fast enough for the integrals to exist. We have 
thus really proven: 

THEOREM: Given two functionals So^], of a finite set of functions such that / and 
So are well-defined on all functions Xi(t), that / > 0, and such that the following functional 
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integrals exist: 5 

Z(g) = J [Vxi] e 



(S [xi]+gl[xi]) 



k\ 



J [Vx^ I[xi\ k e- So[x ' ] < k < N + 1 



Then: 



N 



\Z{g)-Y,Z k g k \<\Z N+l g N+1 \ (11) 



fc=0 



IV. APPLICATION TO FIELD THEORY 



The theorem formulated in the last section is general enough to include the case of 
regulated quantum field theory. As an example consider a self-interacting scalar field theory 
in d < 4 dimensions defined on a lattice in a box: the functions being integrated over are 
4>x{t), the boundary conditions are chosen for example to be periodic, and the functionals 
are: 

^° xglat. 

r <tt ^ 

ccElat. 

The functional integrals Z and Z k exist, and thus the bound (fTTjl applies. 

Unfortunately, the perturbation series whose coefficients are Z k is not the perturbation 
series we are actually interested in. As we remove the lattice and box regulators it will 
almost certainly diverge, and even if we do not the strong dependence of the higher orders 
on extremely high energy physics violates our sensibility. We must instead introduce renor- 
malized perturbation theory before we can produce physically meaningful results. However 
the renormalized coupling enters the functional integral in a subtle way: we demand that the 



o 



bare couplings Zq 6 , mo, and go depend on both the regulator and the renormalizec 
and mass in such a manner as to satisfy Wilson's renormalization group equations 



coupling 



10|. This 



functional dependence can be extremely complex, and there is absolutely no guarantee that 



5 Here recall from footnote 2 that I mean they exist after being discretized. This theorem will be agnostic 
about the limit as the regulator is removed, but for Gaussian So the limit can usually be taken and the 
theorem will be then true for the limits. 

6 The field renormalization coefficient, not the Oth order of perturbation theory! 
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the simple manipulations of the previous section will generalize. Moreover for renormal- 
ized perturbation theory even the standard large-order arguments I sketched in section II 
run into difficulties related to "renormalon" singularities that can dominate the instanton 
contributions to large orders 



Nonetheless I expect that the situation is better than it appears. The problem with our 
calculation is that we are attempting to use a single renormalization scheme to produce a 
perturbation series for all physical processes. The lesson of the renormalization group is 
exactly that this is how we should NOT proceed. Instead we should consider the functional 
integral expression for a specific "infrared-safe" observable at a specific energy scale, and 
then choose a renormalization scheme that is appropriate for that scale. It is only this 
perturbative expansion that really must obey the error bound, and it seems plausible to me 
that it will. I believe that by carefully examining this perturbation series and its relation to 
the bare perturbation series, it should be possible to generalize my argument above to the 
renormalized series, but I have not seen a simple way to do so. 

V. ENERGY LEVELS AND CONNECTED S-MATRIX ELEMENTS 

So far I have discussed Euclidian transition amplitudes at finite time, but in fact these 
are not the objects that are normally studied using perturbation theory: indeed they have 
a "large time" infrared issue as I will now argue. 7 For most experiments the relevant time 
interval in the transition amplitude is much larger than the microphysical timescale of the 
system in question. Moreover this time-dependence becomes stronger and stronger at higher 
orders of perturbation theory, essentially because at higher orders we can have more and 
more disconnected Feynman diagrams and each diagram will be proportional to a power of 
T with the power given by the number of connected components of the diagram. This means 
that regardless of the smallness of g, for large T the perturbation series will be useless. This 
is consistent with the theorem proven above since the bound becomes trivial in the large T 
limit, but it also means that the theorem is not quite what we would want. Thus even in 
the case of simple quantum mechanics we need to do some sort of infrared resummation to 
get a useful result. 

7 This problem was emphasized to me by E. Witten. 
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The observables that are usually approximated by perturbation series are energy levels 
and connected S-matrix elements, so it is for these that we would really like to find a non- 



trivial bound. A colleague and I have spent some time thinking about this, [12J but so far 
we have no decisive result. I will say a little about how the problem can be approached, and 
quote an intermediate result we have obtained. 

The infrared resummation is easy to perform in these cases; it amounts to simply taking 
the logarithm of the perturbation series. For example if Z = tr [e~ TH ] and Zq = tr [e _Tf/ °] , 
then to find the shift in the ground state energy (assuming a discrete spectrum in the vicinity 
of the ground state) we have: 

E — E — — lim ^log|- (12) 

A similar result holds for the generating functional of connected correlation functions. 8 
We would thus like to find a way to bound the error of the perturbative expansion for the 
logarithm of a path integral. The result that we might expect to be true is that the previous 
result continues to hold but with with only connected diagrams appearing on the right-hand 
side of equation [TTJ Indeed using methods similar to those of section three (and under the 
same assumptions) we have shown that: [l^ 

\log Z(g) - logZo| <g^- 

This result is quite interesting, for example if we combine it with equation [TJ] we see that 
the absolute value of the ground state energy shift is bounded by the absolute value of the 
first order correction. 9 The reader can see that the fraction on the right-hand side cancels 
all disconnected diagrams to first order in perturbation theory, so the large T behaviour of 
both sides is linear and the bound is now nontrivial in the large T limit. Unfortunately with 
no other restrictions this bound does not generalize to higher orders; in fact my colleague 
and I have found (admittedly unphysical) counterexamples. Nonetheless it does seem to be 



(13) 



8 See for example chapter 4 of [13]. Correlation functions are then related to the S-matrix via the usual 
LSZ argument. 

9 The reader might be surprised by this since there are well-known examples where the first correction to 
the ground state energy vanishes, but these always involve either fermions or an interaction unbounded 
from below, so the assumptions of the theorem are violated. In fact this bound for the ground state energy 
also follows from a simple variational argument, as was pointed out to me by L. Susskind. 
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obeyed in a variety of interesting examples, so we expect that with a sufficient restriction it 
may still be true. 



VI. CONCLUSION 

The situation is thus mixed. For short-time transition amplitudes in simple quantum 
mechanical systems we have a very satisfying result; in particular it assures us that as long 
as we find successive terms in perturbation theory to be smaller and smaller, we are getting 
closer and closer to the exact answer and should keep computing higher orders. Given the 
potentially divergent nature of the series, this was not at all obvious beforehand. Moreover 
we see that the smaller the coupling, the more accuracy we can achieve before the series 
begins to diverge. This justifies assigning the label of "size of non-perturbative effects" to 
the smallest term in the series. In fact the theorem has as an obvious corollary the usual 
statement that the perturbation series is asymptotic in the mathematical sense, so in this 
case we have managed to derive both parts of the lore stated at the end of section II. However 
so far there is no satisfactory generalization to energy levels and connected S-matrices, as 
well as no obvious way to take into account renormalization. For the first of these I have 
quoted a promising initial result, but it is fair to say that both questions remains open. 
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APPENDIX A: PROOF OF THE IDENTITY 

In the text we used the following identity: 




(Al) 
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To derive it, observe: 10 

ATlg-V - (-x) k = iV' - V [ ~ X) 

k=0 v ' v ' k=N+l 

N , 



X 



AT! e ~ x - V 

(-x) N + 1 [ ^ k\ 

N\ ( .A (-x) 



\ k=0 



N '-xY^ 



= n\ e x V 

Moreover using multiple integration by parts, we can show that: 



dtre xt = -e x + '-e x - ... - -— e x + 



xx 2 x 3 "' (-x) N+1 (~x) N+1 

By inspection these two sums are equal and thus we have proven the identity (lAil . 
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